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ABSTRACT 

We present a new method to compute the first crossing distribution in excursion set 
theory for the case of correlated random walks. We use a combination of the path 
integral formalism of Maggiore & Riotto, and the integral equation solution of Zhang 
& Hui, and Benson et al. to find a numerically robust and convenient algorithm to 
derive the first crossing distribution in terms of a perturbative expansion around the 
limit of an uncorrelated random walk. We apply this methodology to the specific case 
of a Gaussian random density field filtered with a Gaussian smoothing function. By 
comparing our solutions to results from Monte Carlo calculations of the first crossing 
distribution we demonstrate that our method accurate for power spectra P{k) cx k" 
for n = 1, becoming less accurate for smaller values of n. It is therefore complementary 
to the method of Musso & Sheth, which will therefore be more useful for standard 
AGDM power spectra. Our approach is quite general, and can be adapted to other 
smoothing functions, and also to non-Gaussian density fields. 
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1 INTRODUCTION 

The growth of large scale structure in our Universe, and the subsequent formation of galaxies, is a challenging, non-linear 
problem. Small perturbations in the matter density in the early Universe grow via gravitational instability and eventually 
collapse to form dark matter halos. Haloes grow through merging with other halos in a hierarchy of ever increasing mass. This 
hierarchy of dark matter halos provides the cosmic web within which galaxy formation takes places. It is therefore crucial to 
understand the formation of dark matter halos in order to accurately model the formation of galaxies ( Baugh|2006 1. 

N-body simulations have developed into an extremely powerful tool to explore the formation and evolution of dark matter 
in the highly non-linear regime. However, they remain computationally expensive, making them impractical for surveying 
dark matter parameter space for example. An alternative approach to computing the formation of dark matter halos was first 



described by Press & Schechter ( 1974 1, and later developed extensively by Bond et al. ( 1991 1 and Lacey & Cole ( |1993) , using 
excursion set theory. This extended Press-Schechter (ePS) methodology provides a highly successful statistical description of 
the formation histories of dark matter halos. For a recent review of excursion set theory see |Zentner| ( |2007l ) . 

The formation of halos has been studied using many different approaches, including numerical approaches such as Monte 



(19911; 



Carlo (e.g. Bond et al. (19911; Paranjape et al. (2012l), and analytical approaches (e.g. Bardeen et al. (19861; Bond et al. 
Maggiore & Riotto ( 2010a[ ); Paranjape & Sheth (2012l). Of these analytic methods, the two most common are peaks 



Bardeen et al. ( 1986 1, and the excursion set theory ^Bond et al. ( 1991 1 



theory. 

The excursion set theory used in the ePS technique follows random walks in the space of overdensity and variance (of the 
density field), an associates halo collapse with the first-crossing of a barrier by a random walk. For the case of a constant, or 
linear barrier and for random walks in which each step is uncorrelated with previous steps, an analytic solution can be found 
for the distribution of variances (and, therefore, masses) at which halo collapse occurs. This allows the halo mass function to 
be found directly for the simple case of a Gaussian-random density field and a variance derived using a filter that is sharp 
in fc-space. Numerical methods exist to solve the problem numerically in the case of non-linear barriers ( Zhang fc Hui||2006[ 



^ Recently, jParanjape et ah|2012^ and jParanjape fc Sheth] ( |2012^ combined the excursion set and peak theory methods to construct an 
excursion set theory for peaks. 
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Benson et al.|2012 1, but non-Gaussian density fields and/or different filters, botli of wtiicti cause successive steps in ttie random 
walk to become correlated, are more problematic. 

Many attempts have been made to use path integral techniques- 



-a powerful tool for studying stochastic and random 
processes — to solve problems in excursion set theory (e.g. Maggiore fc Riottol (|2010a|b|^ 



Ma et al. (20111; Adshead et al. 



(20121). In this work, we modify the path integral method of Maggiore & Riotto (2010a I and combine it with our modified 
integral technique for solving the first crossing distribution, which is based upon the work of Zhang & Hui ( 2006 1 and which 
we introduced in our previous work, [Benson et al.| ( |2012| ). Through this combination, we develop a new approach to solve the 
first crossing problem in more general cases, and specifically address the case of a Gaussian filter. 

The remainder of this work is organized as follows. In §2.2| we review how the path integral approach is used to find the 
probability distributions of trajectories and describe how to add a perturbation term representing the effects of the filter. In 
i|3] we apply these methods to find the first crossing probability, while in !|4] we demonstrate the method using a Gaussian 
filter. Finally in ^we discuss our results. 

We also include two appendixes. Appendix |X] presents a detailed calculation for the second-order contribution to the 
perturbation term and outlines how this same procedure can be used to derive higher order perturbation terms. Appendix [B] 
presents our numerical technique for solving the first crossing distribution which we have found to be robust against numerical 
error. 



2 PROBABILITY DISTRIBUTION OF TRAJECTORY 



2.1 Introduction 



As shown by|Maggiore fc Riotto' f2010aVMaggiore fc Riotto|(|2010"b|, and'Maggiore fc Riotto] (|2010c[) (hereafter MRIO) in their 



series of works, path integral methods provide a powerful tool for studying the formation of dark matter halos in the context 
of excursion set theory. In their works they used this method to address the issues of non-Markovian trajectories, moving 
barriers, and non-Gaussianity in excursion set theory. Additionally, Zhang & Hui ( 2006') introduced an elegant method for 



solving the excursion set theory for moving barriers given knowledge of the probability distribution function for overdensities 
as a function of variance. 

All of these methods have limitations, which we will discuss further in !|5] In this work we modify the methods of MRIO 
and Zhang & Hui ( 2006 1 and combine them to solve the excursion set problem in the general case with greater accuracy. 



Specifically, our approach differs from that of MRIO as we find the probability distribution function for trajectories without 
consideration of the barrier, and then combine this with the method of [Zhang fc Hui| ( |2006[ ) to find the first crossing distribution 
for an arbitrary barrier. 

In this section we review how the probability function for a trajectory can be found utilizing path integral techniques, 
and how one can find the relevant first-order perturbation term for different filters. In Appendix [A| we present a calculation 
of the second-order perturbation term, which, for many cases, will be sufficient. 



2.2 Calculation of the Path Integral 

We begin by applying the path integral method to the case of a Gaussian random density field with a sharp fc-space fitlering 
and demonstrate that in this case the well-known analytic solution is recovered. 

We consider the density perturbation, 5(x), at point x in the universe, defined as 



= pM_I, (1) 

In ePS theory we are interested in a smoothed density field, and so we smooth this density field with some filter function. 
Assuming a Gaussian-random density field, if the field is smoothed with a sharp fc-space filter (i.e. a filter that is constant 
below some critical wavenumber, and zero otherwise) the resulting random walk in 5 as a function of smoothing scale will be 
a Markovian process. This has a well-known analytic solution for the case of a linear barrier 
2002 1. 

After smoothing, the variance of density field is 



Sheth[p98 Sheth fc Tormen 



S = J dink A^ik) \W\k,R)\, (2) 

where W{k,R) is the filter function in Fourier space, R is the smoothing scale, and A{k) is the dimensionless power 
spectrum. In the path integral approach, we assume that S{S) evolves stochastically with "time-like" S. By definition, {S(S)) — 
0, and for Gaussian fields the only non-vanishing connected correlator is the two-point correlator, {S{Si)5{S2)) ■ 
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We wish to find the probabihty distribution of trajectories which start from S — So, 5 — So and finish a,t S = Sn, S = 5n. 
We imagine discretizing our domain, [So, Sn] into n equal sized parts. Each step is therefore e = AS/n. The trajectory travels 
through S and in each step Si takes on a specific value. We will integral over all Si's from —oo to oo to find the total probability 
for the trajectory given its starting and ending points. In the limit n — >■ cxa we will recover the continuum solution. We define 
the probability density in the space of trajectories as 



W{So; Si,...,S„-S„) = (SniSiSo) - <5o) ■ ■ ■ <5b(<5(S„) - 5„)), 



(3) 



where, for clarity, S is the density fluctuation field and Sn is the Dirac delta function. The probability distribution is 
found by integrating over all possible paths and has the form 



P{So', So; Si, ... , Sn~i; S„; Sn) — 



dSi 



dSn 



X W{So;Si,S2,S3,. . . ,S„-i;Sn). 



(4) 



For the case of a Gaussian random density field and a sharp fe-space filter this probability distribution depends only on 
the starting and ending points, as we will show below. The function W{So; Si, ... , S„-i; Sn) can be expressed in terms of the 
connected correlators of the theory. 



W{So; Si,..., Sn-i; S„) = 
or more practically. 







J — oo 


J —oo 



27V 



2tt 



(5) 



W{So; Si,. .. , S„-i;Sn) 



I 1 = 1 p=2 ^' 11=1 ip = l 



27r 



2tv 



(6) 



Here {Si^ ■ ■ ■ Si^) denotes the connected n-point correlator. For a Gaussian random field density it is possible to show that 
the only non-vanishing connected correlator is the two-point correlator, {S{Si)S{S2)) . Then for the Gaussian case we find. 



W{So;Si 



, Sn-i;S„) 



dXo 
2tv 



dXn 

2tv 



f n ^ n 

exp < i^^XiSi — - ^ XiXj{SiSj) 



(7) 



2.3 Gaussian Random Field with Sharp fc-Space Filtering 

It is well-known that for a Gaussian-random density field and sharp fc-space filter the probability density function for trajec- 
tories is 



P{So,So = 0;Sn,Sn) = 



: exp 



{Sn - So)^ 



y/2TTSn I 2Sn 

Using eqns. Q and Q we should therefore be able to reproduce this known result. We find 



(8) 



P{So, So; Sn, Sn) 



dSi 



dSn~l 



(9) 



^ . . . ^ exp i i ^ \S, >^^^j(S,Sj) 

{ 1 = 1 i,i = l 

where (SiSj) is the two point correlation function. MRIO showed that for a Gaussian random field and sharp fc-space filter 



the two-point correlator has the form. 



{S{S,)S{Sj)) = min(S„5,) = emin(j,i). 
To solve the integrals in eqn. (|9| we will use the following result (see, for example, Altland & Simons poio| . 



(10) 



(11) 



We first integrate over the A's in eqn. Q by setting bi = Si, Aij = (SiSj), and Xi = Xi to giv^ 



^ For a detailed calculation, see MRIO. 
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d<5i...y d5„_i^^— p^exp|--^(5,+i-50'|- (12) 

We now employ eqn. ||9| again, this time with Xi — Si. In this case, the elements of matrix A are Ai,i+i = = — 1/e 

and Ai^i — 2/e, and the determinant of matrix A is ne^~". The elements of vector b are: bi = 6„-i = — e^^ and the other 



elements are all zero. Using eqn. 1 11 1 to solve eqn ( 12 1 we find 



P(5o;So = 0;5„;5„) = — ^exp|-^^^V^V (13) 
V27rne I 2ne J 

Using the fact that ne = S'n this reduces so, 

P(5o;So = 0;5„;S„) = ^^exp|-^^^^^^|, (14) 
which is independent of our discretization and agrees with the well-known analytic solution. 

In the above we have assumed So ~ 0, but in general, one may be interested to have the probability distribution of 
random walks starting from arbitrary Sq. The two-point correlator is unchanged in this case, and so the above calculation 
holds, but now ne = Sn — So, and so 

For solving the excursion set problem it will be useful to calculate the probability of a trajectory which starts at point 
5*0 = 0, 5 = and then passes through the point S — Si, S — Si before ending at S — Sn, S — S„. We will use this result in 
section |3] In this case, we simply do not integrate over Si in eqn. ([£]), 

/-j-OO /' + 00 /' + 00 /' + 00 

dSi... dS^-i d(5.+i.../ dSn-iW'"'\So;Si,--- ,Si,--- ,Sn;S„). (16) 
-oo J — oo J — oo J — oo 

It is then straightforward to show that 

P^A A Q A Q ^ 1 / {S^ - Sof \ ( (Sn - S,f \ 



2.4 Gaussian Random Field with Arbitrary Filter 



In this section we use the prescription introduced in §2.3| to calculate the probability function for trajectories with an arbitrary 
filter. For an arbitrary filter it is not possible to solve the necessary integrals analytically. Therefore, we adopt a perturbation 
theory approach to this problem. In this section we find the first-order solution, while Appendix |A] gives the second-order 
term. Finding higher order terms is possible, but in many cases second order will be sufficient. 

Without loss of generality will will assume that So = and So =0. As we are still considering a Gaussian random field 
all correlations except the two-point correlation will vanish. We therefore assume that a small perturbation, A{Si, Sj) = Aij, 
is added to the two-point correlation. We refer to Aij as the "filter effect" term. We then have. 



P{5o;Sn;Sn) 



dSi . 



dSn- 



+°° dAo dA^ 
2n "' 2n 



exp < i ^ XiSi — - ^ XiXj (min {Si,Sj) + Aij 



(18) 

Calculation of Aij for a given filter is discussed in § 2.5 Defining di = d/dSi, we can rewrite equation 18 retaining only 
first order corrections, as 



P{So; Sn, Sn 



/+00 /'+00 
d5i... / d5„„iVK^'"(5o;5i,--- ,Sn;Sn) 
-00 J — 00 



dSi 

+ 2 E ^^^1 



dSi 



dSn-id,d,W'^{So;Si,--- ,Sn;S„), 



(19) 
(20) 



where. 
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" J ^ ■ ■ ■ ^ """"P I* ^ ~ ^ ^ ^"^^ J 



(21) 



Note that Aij — Aji and A^^ = 0. Since the two point correlation lunction is symmetric and we always have (SiSi) — 
mm{Si, Si) — Si, the perturbation to the two-point correlation function should be zero for Si = Sj. Using An — and the 
symmetry of the perturbation term one can write, 

^ n ^ n — 1 rt — 2 n — l 

- J2 ^^Adj = 2 E + ^ ^ ^^^^ 

— 1 i — l i — l J — i + 1 

MRIO argued that YJ^Zi ^indid„ gives an integral over an intermediate time, say Sj, and X^ILi^ SjLi+i produces a 
double integral over intermediate time, 5*. They therefore refer to these respectively as the "memory term" and "memory of 
memory term" . We will adopt the same labels in the remainder of this work. 

Thus, 

P{So; Sn;Su) = P"'"(5o; 5„; S„) + P"^""(5o; 5„; S„) + p"— -°-(5o; 5„; Sn). (23) 
The probability distribution P^™ was derived in j]2.3| We now consider the first perturbation term, p™'=™: 



P'""'"(5o;5n;S,0 = EA,„a„ / dS^ . . . d5„_iftH/«^'"(5o; 5i, ■ • • , 5™; S„). (24) 



In the continuum limit eqn. ( |24[ | this becomes an integral over time, S. We proceed by working with the discretized form 
and at the end will take the limit of n — >■ oo to find the continuum result: 

/' + 00 /' + 00 /' + 00 /' + 00 

pmcm^^^. g^^-j ^ ^ Aindn Hm / d5i . . . d(5,;_l / dS^+l ...I d(5n-l 

^^^'^^ J ~oc J — oo -/ — oc J ~oo 

X WS'"(5o;5i,--- = ,<5„;S'n) 

= ^ Aindn lim / d5i . . . d(5,;_i / d^^ + i . . . / dJn-i 

Ai„9„ lim / dSi . . . dSi-i / d5i+i . . . / d(5„_i 

f-f ^c-^ + ^o J-oo J-oo i-oo i-oo 

X W^'^'"(5o;<5i,--- ,5c; SO x W''\5o■,S^+l,■ ■ ■ ,<5n; (S„ - SO), (25) 



1 

2^ 

which can be simplified to give 



n-l 

P'"""(<5o; 5„; S„) = V A.„a„ lim P'^'"(5o; = 5,; SOP'^'^l^. = S,; S„; (S„ - SO). (26) 

i—l 

It is obvious that for i = 1, - • ■ ,n — 1, 



dn lim 

5c— > + o 



dSi... dS.-i dS.+i... d5„_i xH/'^'"(5o;<5i,--- ,c5« = 5c,--- ,5n;S„) = 0. (27) 

-oc J — oo J — oo J — oo 

Therefore, P"""^{So; (5„; S„) = 0. In the same way, one can prove that pmcm-mom^^^. g^. ^ implying that the whole 
perturbation term becomes zero as 5c — >■ +oo. Consequently, it is possible to write the solution, when no barrier is considered, 
for all filters in terms of min (Si, Sj) + Aij, and all perturbation terms to the probability function become zero. 

Now we again assume that the trajectory passes through point S — Sk, S = Sk- We will show that in this case the choice 
of filter does matter. For the memory term we find 



P'""'"(5o;<5n;S„) = Afe„9„ / dS^... dS^.^ d^+i • • • / d5„_i9feVK^"^(5o; 5i, • ■ • S„), (28) 

J — oo J — oo J — oo J —oo 

where Akn is the perturbation term and dk = d/dSk and 5o = 0. This gives, 
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Figure 1. The FE for a Gaussian filter as given by eqn. i32t. Blue, red, and green lines show the FE for Sj = 1, 3, 5 respectively. 



P"""^(5o; S„; S„) =Aknd„dk 



dSk-i / dSk+i 



dSn 



1 ( 1 



-Akndnd) 



— SlSn — SkiSn — S„ + Sk) + SkSn{Sn + Sk) ( {Sk ^ Sq)^ \ ( {Sn — Sk)'^ 

= - ^fcn ^-3 ,n . _ z-ttt: X exp <^ — } X exp •' 



27rS^/'(S„-Sfe)5/2 



25-,. 



2(5„ — Sk) 



(29) 



For the memory-of-memory term we have, 



-imem — mem 



{5o;5„-Sn) = Aij^ / dSi 



dSk-i / dSk+i . 



d5„-ididjW^'^{5o;Si 



1 \ Sn 



(30) 



Therefore, all terms of p™°™ vanish because there is at least one derivative with respect to Si for which i ^ k which 
cancels the integral over Si. Taking the limit Si — >■ +<x leaves us with zero. 



2.5 Calculating the Filter Effect 

We have not yet discussed how to calculate the filter effect (hereafter FE), A{Si,Sj) = Aij. A particular advantage of our 
method compared to its ancestors is that the terms that contribute to the FE are completely independent of the perturbation 
terms. This makes it very easy to develop solutions for alternative filters which can be used in eqn. ( |29[ ). 
The non-perturbative part of the two point correlation function was min (5*;, Sj). Therefore, 



{S.S,) 



(31) 



One can use the above formula to find the perturbation formula for any filter of interest. Note that, for most cases there 
will be no analytical solution, and the filter effect must be computed immerically for each point. 

To be specific, we will now calculate the FE for a Gaussian filter. For simplicity, we consider a power spectrum of the 



form P{k) oc k" with n — 1. Using eqn. (481 and assuming i < j we find 



3S.S, ~ - 2^/{SfS~) 

Si -\- Sj -\- 2^ SiSj 



(32) 



where GF stands for "Gaussian filter" . Note that A*]^ in eqn. ( 32 1 is not symmetric in its arguments as we have explicitly 
assumed that i < j. 

In Fig.[l]we show the FE at three points. As expected, for small 5* the FE is less than min(S'i, 5*^) and can be treated as 
a perturbation. For large variances the FE increases, so for large S the perturbation calculation will begin to break down. 
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3 FIRST CROSSING PROBABILITY 



Zhang & Hui ( 2006 1 introduced an elegant method to find the first crossing probabihty in excursion set theory for arbitrary 



barriers. In our previous work ( [Benson et aL|[2012[ ) we modified this method to obtain more accurate answers for moving 
barriers. In this section we will incorporate the filter effect into this integral formulation. The method of Zhang & Hui (20061 



begins by considering the integral form of the first crossing problem, which simply states that the total probability at any 
variance 5* is the sum of the probability that a trajectory has crossed the barrier at some smaller S, and the probability that 
the trajectory is currently below the barrier having never crossed at any smaller S: 



f{S')dS' 



i3(S„) 



PB{S,Sn)d5 



(33) 



where Pb{S, Sn) is the distribution of 5 at Sn accounting for the absorbing barrer and is given by 



PB{Su,S„)^P{5n,Sn) 



G(S„,<5„,5")d5'. 



(34) 



In this formulation, eqns. ( 33 1 and ( 34 1, G{S„, Sn, S')dS' is the probability for a trajectory which crosses the barrier for the 



first time between S' and S' + dS' and approaches the point {S„, S„)- It is calculated by multiplying the fraction of trajectories 
which cross the barrier between S' and S' + dS' for the first time, f{S')dS', by the probability, P{5o; So; S' , B{S'); Sn', Sn), 
of a trajectory which starts from point {So =0,5 = 0), crosses the barrier at point (S = S',5 = B{S')), and finally ends at 
point {S = Sn, S = 5n)- Here, B{S), is the moving barrier. In general one can calculate the second term of equation 34 using. 



nfQ X Q'^ — ffc'\ ^(^oi ^ , B{S );S„; Sn 



(35) 



where P{So; So; S' , B{S')) is probability of trajectory which starts from point {S = 5*0,5 = So) and ends at point 



(5* = S',S = B{S')). In the case of a Gaussian-random field and a sharp fc-space filter one can divide eqn. (171 by eqn. (14 1 
and obtain 



G{S,S, S',B(S')) = fiS')- 



: exp 



VMS-S') 

For the case of a Gaussian-random field and a generic filter one can use eqn. ( 29 1 to obtairj^ 



G{S,S,S',BiS'))=f{S') 



V27r(5 - S' 



: exp 



^ -B(g')5 - S'jBjS) -S + S') + BiS')B{S){S + S') J iB{S) - B(S'))' 



^27rS"2(5-5''; 



2{S-S') 



(36) 



(37) 



We define the second term in the above equation to be f(S')A{S', S)'y{S, S, S' , B{S')). As we have shown in Benson et al 



(20121, to discretize this equation and find the first crossing probability, f{S), the next step is to find the following integral. 



$(5, S') = A(5', S)7(S, S, S', BiS'))dS = -A(S, S') x M|1|_£MJ ^ ^ ^ . (33) 

J —00 



y/2n S'^^iS - S')'-^ 



2{S-S') 



We can now discretize eqn. (33 1 and solve for f{S). In our previous work we showed in detail how one can discretize this 
equation. Briefly, we have 



/(«') 



erf 



B{ S)-B{S' ) ] , 
V^S-S') J 



dS'. 



(39) 



We assume that AS = Sma^/N. For eqn. (|33[) we find 



^ Note that, in the general case, one should consider not the fraction of trajectories passing through point {B[S], S) which end at point 
{5n,S„), but the fraction of trajectories which pass through point (SIS'], 5) having never crossed the barrier for smaller S which end 
at point {Sn,Sn)- For a Gaussian random field with a sharp fc-space filter, these two are equivalent, as the path of any given trajectory 
for 5' > 5 is independent of its path for 5" < S. For other filters, there is some correlation with the path for S" < S. Therefore, while 
eqn. 13611 is exact, eqn. (|37|l is an approximation, obtained by ignoring the correlated nature of the random walk. 
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E 

j=0 



Since /(O) = we obtain, 



AS 
2 ■ 



(40) 



1 - erf 



B{S,)-B{S,) 



AS' 



1 - erf 



Bis, 



1 - erf 



BiS,)-B{Si 



2$(5,,S0 I AS. 



a/2(S, - SO 

For all barriers of interest erf|B(Sj) - B{Sj)/^2{Sj - Sj)| = and $(5^, Sj) = 0, such that 



1 - erf. 



- J2 fiS^) 1 - erf - 2$(S„ so AS 



\/2(S, - SO 



2 

AS' 



In the above, $(Sj,SO plays the role of a filter correction term, and is calculated using eqn. (38 1 



(41) 



(42) 



Here we have computed up to first-order perturbations. In Appendix |A] we derive results up to second-order which is 
sufficient for most purposes (although this will depend on the choice of filter). We address numerical issues in solving for f{Sj) 
in Appendix [B] 



3.1 Gaussian Random Field with Gaussian filter 

So far, our calculations are independent of the choice of filter. In this section we compute the first crossing distribution for 
a Gaussian filter applied to a Gaussian-random field. For a Gaussian-random field both the unconditional and conditional 
probabilities should also be Gaussian, independent of the choice of filter. We use this case to illustrate that the path integral 
formalism gives the expected answer. Of course, the true power of the path integral formalism will become relevant when 
applied to non-Gaussian fields. The Gaussian filter is convenient as its two point correlation function can be calculated 
analytically and, unlike the sharp fc-space filter, its smoothing region in configuration space is well defined and localized. In 
this case, the variance, S, is defined as 

where 



S = {S\R)) = 1^ P{k) W\k, R) dk, (43) 



Wik,R) = e^p{~^} (44) 



is the Fourier transform of our filter which in real space has the form 

(27ri?2)3/2 '''^^{~2R^ 



W{r,R)= ,^^L:,,, expj--^^. (45) 



For a scale-free power spectrum, P{k) — Ak", 



S ^ {S\R)) ^ (^^^ R-^-" , (46) 

and two-point correlation function is 

<~))^j^.r(i±^)(| + |)"^, 

where here F is the usual Gamma function. For n — 1 therefore, we find the following relation between the variance and 
the two-point correlation function: 
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(S.S,)^ (48) 

5*1 + Sj S-iSj 

For general exponent n this becomes: 

{5.5,) = I^J11±^]^ . (49) 

Setting i = j gives us the variance of point Si, (SiSj) — Si, consistent with that found for a sharp fc-space filter for which 
= min(Si, Sj). 



4 RESULTS 

We can now use our results to compute the first crossing distribution for a Gaussian filter. MRIO also employed the path 
integral method account for the filter effect as perturbation term. They used an approximate form for the two-point correlation 
function for both a top-hat filter in real space and a Gaussian filter. Specifically, they proposed the following function: 



1 — K Sc 



- exp 



Si 
2S 



+ 



2V2^ sy- 



-r 



51 
25 



(50) 



/27 

where r{0,z) is the incomplete Gamma function, Sc ~ 1.686 is the barrier for spherical collapse. MRIO computed the 
value of K for top-hat and Gaussian filters, and a ACDM power spectrum, finding that k ~ 0.4592 and 0.35 respectively with 
a weak dependence on scale. Because MRIO computed k for a ACDM power spectrum it is not directly comparable to the 
power-law cases we consider here. Furthermore, for case of power-law, we have compared the form of our expression for Aij 
and find that it is not well approximated by MRlO's fitting function for any value of k. Therefore, we instead compare our 
solutions to the results of Monte Carlo simulations of the first crossing problem. 



Peacock & Heavens ( 1990 1 used an alternative method to find the first crossing probability distribution for correlated 



random walks, proposing an approximate formula for the Gaussian filter. Later Paranjape & Sheth (20121 proposed another 



approximate formula and argued that their formula worked well for a large range of S. More recently, MS12 used this method 



to propose an improved approximation, which is similar to that of Peacock & Heavens (19901, for generic filters as well as 
moving barriers. 

Fig. [2] compares results from this work with those of MS12, and a Monte Carlo calculation of trajectorie^ We use 
Robertson et al. (2009l's prescription (see also Bond et al.||1991 l for Monte Carlo simulation of random walks. Results are 
shown for two different values of the power spectrum index, n. For n = -1-1.0 (left panel) our results are in very good agreement 
with the Monte Carlo results for v > 1 when we include up to the second-order perturbation term. For smaller 1/, our result 
begins to diverge from the Monte Carlo result, albeit slowly. Higher order perturbation terms could be included of course, 
which would extend the validity of our result to smaller i/. The approximation of MS12 performs poorly for n = -f 1.0, differing 
significantly from the Monte Carlo solution for h' < 6. For n = —1.2 the situation is reversed. Here, the approximation of 
MS12 works extremely well (as MS12 showed), while our approximation disagrees substantially with the Monte Carlo result. 

Our method applies equally well to non-constant (a.k.a. "moving") barriers. Figure [s] compares results from this work 
with those of MS12, and a Monte Carlo calculation of trajectories for a linear moving barrier with the form of. 



B{S) = Sc{l + ^). 



(51) 



The figure shows the result for two choices of a, +0.5 (dashed lines) and —1.0 (solid lines), and two different power 
spectrum, n — +1.0 (left panel) and n = —1.2 (right panel). For both power spectra, the approximation of MR12 and this 
are in close agreement for a < (although MS12 remains the better match to the Monte Carlo results for n — —1.2), but as 
a is increased they diverge from each other, such that for n = +1.0 (left panel) our results are in very good agreement with 
the Monte Carlo results, while for n = —1.2 the approximation of MS12 works better (as MS12 showed). 



We cannot compare with the results of MRIO in this case. Their mass function is expressed in terms of a parameter, k, which 
parameterizes the perturbation, Ai^-. We find that the parameterization used by MRIO — derived for a ACDM power spectrum — is not 
a good match to the perturbation derived in this work for any value of re. Note that MRIO used transfer function for calculating re but 
in this work just the power-law is considered. For similar reasons, the work of |Corasaniti fc Achitouv| ( ,2011^ and |Achitouv & Corasaniti| 
1 2012 I are also not comparable to our current examples. 
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Figure 2. The dimensionless halo mass function vf{v) as computed using different approximations, for a Gaussian-random density field, 
smoothed with a Gaussian filter, and for two different power spectral indices: n = +1.0 (left panel) and n = —1.2 (right panel). The dark 
blue line shows the results from a Monte Carlo calculation, while red and magenta lines show results from this work including up to the 
first and second-order perturbations respectively. Black lines show the approximation of MS12. For reference, the green line shows the 
exact result for a sharp fc-spacc filter. 



n = l n = -1.2 
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Figure 3. The dimensionless halo mass function yf{v) as computed using different approximations, for a Gaussian-random density field, 
smoothed with a Gaussian filter, and for two different power spectral indices: n = +1.0 (left panel) and n = —1.2 (right panel), for the 
case of a linear moving barrier of the form B{S) = 5c(l + oiS/S^). Dashed lines show results for a = +0.5 while solid lines show results 
for a = —1.0. The dark blue line shows the results from a Monte Carlo calculation, while red and magenta lines show results from this 
work including up to the first and second-order perturbations respectively. Black lines show the approximation of MS12. For reference, 
the green line shows the exact result for a sharp fe-space filter. 



5 CONCLUSION AND DISCUSSION 



Many attempts have been made to account for correlated random walks in excursion set theory, whether induced by the filter 



effect or by intrinsic correlations in the density field. For example, Peacock & Heavens ( 1990 1 computed approximate solutions 



for correlated random walks, specifically for a Gaussian filter function. This approach was later improved (Paranjape & Sheth 



2012[|Paranjape et al.|2012[[Musso &i Sheth|2012||Musso et al.|2012[ ) to work with moving barriers. While direct Monte Carlo 
simulation can account for these correlations it is inefficient and therefore impractical for application within semi-analytical 
models which use ePS theory to construct merging histories for dark matter halos. MRIO employed a path integral method 
to include the filter effect as a perturbation term. They used an approximate function for the two-point correlation function 
arising from both real-space top-hat, and Gaussian filters. In their later works they have proposed a series solution for moving 
barriers with sharp fc-space filters. It is difficult to include both the filter effect and a moving barrier in the approach of MRIO. 
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Later works ( Corasaniti fc Achitouv|2011 1 improved MRlO's method by considering a linear barrier in the perturbative path 
integral method, but still are not applicable to an arbitrary moving barrier. 

The solution for arbitrary moving barriers with uncorrelated random walks is well understood and has been studied 
for many years ( Zentner]|2007[|. In particular, the integral method of Zhang & Hui (20061, and its numerically more robust 
modification ( Benson et al.|2012 1 , work well in this case. In this work we have combined the power of the path-integral approach 
to excursion set theory with the integral-method solution for arbitrary barriers to develop a perturbative solution to the first 
crossing problem applicable to arbitrary filters and barriers. The perturbation terms can be computed directly given the 
functional form of the filter. Our approach utilizes the path integral formalism of MRIO, but uses it to find the unconditional 
probability Po{5, S) — the probability that a random walk reaches the point 5 at variance S — and the conditional probability 
P{Sc, Sc; 5, S) — the probability that a random walk cross the barrier at point S then reaches the point S at variance S — 
considering all possible walks and weighting them equally when constructing the ensemble average. Note that we specifically 
do not use the path integral formalism to compute the probability P{S, S) in the presence of the absorbing barrier. In this 
formalism, the conditional probability is responsible for accounting for correlations between steps in the random walk. Given 
these unconditional and conditional probabilities, the barrier is added to the problem and used to compute the first crossing 
distribution using the method of Zhang & Hui ( 2006 1 . We find that this approach is accurate when compared to the results 
of Monte Carlo simulations over a wide range of variance. 

Our method complements that of MS12, which is designed to be accurate in the limit of smooth random walks (corre- 
sponding to more negative values of the spectral index n). Our approach expands around the limit of uncorrelated steps in 
the random walk (corresponding to "jagged" random walks, and more positive values of the spectral index n). This can be 
clearly seen in Figure [3] in which we show results for n — —1.2 and n = -1-1.0. The MS 12 approach is therefore more useful for 
standard ACDM power spectra (for which the local exponent is negative on scales of interest for dark matter halo formation). 
Finally, we compared our solution with that of MS12 for a linear moving barrier with the form B{S) = 3^(1 + aS/5l) for two 
choice of Q, a = -1-0.5 and a — —1.0. As can be seen in Figure[2j for a < both methods are in approximate agreement with 
Monte Carlo calculations. 

Our method has the advantage of being very flexible. Currently, we have used these methods to derive the unconditional 
and conditional probability distributions (required to solve the first crossing problem) for the simple case of a Gaussian random 
field — demonstrating that this approach correctly recovers the expected Gaussian distributions in this case. In a future work, 
we intend to exploit this flexibility of our method to examine models in which the density field is non-Gaussian. Additionally, 
we note that our approach currently gives equal weight to each possible random walk when constructing the ensemble average. 
Dark matter halos may, in fact, correspond to a special subset of all walks, e.g. those corresponding to peaks in the density 
field ( Paranjape fc Sheth|2012 Paranjape et al.|2012 |, which would require different weights to be assigned to random walks 
in the ensemble averaging. We plan to revisit this issue. 
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APPENDIX A: CALCULATION OF THE PERTURBATION 



In § |2.4| we described in detail the calculation of the first-order perturbation arising from the filter effect. Below we outline the 
calculation of the second-order perturbation term. The same approach can be used to derive yet higher order perturbations 
if necessary. We have found that, for large masses, the first perturbation term is sufficient. Adding the second or higher order 
terms extends the range of validity of our approximations to larger S (smaller mass). 

Again we assume that 5o = and So = 0, and assume a density field with Gaussian-random fluctuations. Therefore, the 
probability of finding the trajectory at {5n,Sn) is, 



P + QC 


r+oc 


/' + 00 


r + oc 


/ dSi.. 








J —oc 


J — oc 


J — oc 


J — oc 



dXo dA„ 
^ 27 



f " 1 " V 

exp < i^^XiSi — - ^ XiXj (min {Si, Sj) + Aij) > 



(Al) 

where Aij is the filter effect(FE) term. By defining di = d/dSi, we can write equation Al including perturbations up to 
second-order as, 



P{So;S„,S„) 



dSi 



d5„_iVK^-"(5o;5i 



+ 2 ^ d5i . . . / 

+ ^ ^ I ^km^^j / d5l . . . / 



dSn-ldSjW^'^{5o;5i,--- ,5-a\Sn) 



(A2) 



where, 



/+00 r-\ 
... / 
- oo J —c 



dAo dAn 



2-K 



2-K 



f n 1 

exp I i ^ XiSi - 7^ ^ ^i^j (min {St, Sj)) 



(A3) 



Note again that the perturbation term has the following properties: Aij = Aji, and Aii = 0. Because the two point 
correlation function is symmetric and {5i5i) = min (Si, 5'i) = Si, the perturbation part of two point correlation function is 
always zero. Using Aii = and the symmetry of perturbation term one can write. 



n— 2 n— 1 



(A4) 



i = l i=! + l 



and. 



- ^ij^kmdidjdkdm ^ I - ^ AinftSn + ^ ^ Aij9i9jj X (-^Afe„9fea„ + ^ ^ Akmdkdm \ . (A5) 

i,j,fc,m — 1 \ i — 1 i — 1 j — i + 1 / \ fc — 1 fc — 1 m — fc + l 

In § |2.4| we demonstrated that after integration all terms vanish when finding the probability of a trajectory that ends 
at {S„,Sn), but when considering a trajectory that passes through {Sk,Sk) as an intermediate step the term with Akndndk is 
non-vanishing for the first-order perturbation. We apply this same argument for the second-order perturbation term and find 
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that the only non- vanishing term is that with A|„9^9^. For higher-order perturbations the same argument can be used to 
show that the only non- vanishing term in the A'^*''-order perturbation is: A^9^ dj^ . 

We define the A'^"'-order perturbation term for the trajectory that passes {Sk,Sk) and ends at {S„,Sn) as: 



Ppci{So; Sk; Sk; Sn; Sn) =j^ Akndndk 



X ; — TTT Gxp , 

— ^iVgiv 1 

''Nl " " ^AlV^SkiSr^-Sk) 



dSi ... I dSk-i / dSk+i 

J —oo 



d5„_i 



[Sk — So)^ \ ( {Sn — Sk)'^ 

X exp < — - — — — — > X exp 



2Sk 



^{Sn — Sk) 



The factor 1/A'^! in Ppor arises from the Taylor expansion of eqn.(All 



To find ^{S,S'), as defined in eqn. (381 we employ eqn. (361 to integrate with respect to 5 from — oo to B{S), and set 
Sk = B{S') and Sk = S' . For the iV*''-order perturbation we find: 



(A6) 



$^(5,5') = J V2^e^p( — ]p^,,{So;Sk;Sk;5„-Sn) 



dS. 



Specifically, the first and second-order perturbation terms are therefore found to be, 



is, S') = -A(5, S') X exp I -Mfl^^ 1 X - 



2iS-S') J ^2tvS'HS~ S')^'' 



(A7) 



(A8) 



and, 



-B'^{S)S'^ + B^{S)B{S')S'iS' + 2S) - B{S){2S' + S){S'^ + B^{S')S - S' S) + B{S')S{3S'{S' ~ S) + B^{S')S) 



(A9) 



V27r5'4(S - S'[ 



respectively. 



APPENDIX B: NUMERICAL ALGORITHM 



In our previous work( Benson et al. 20121 and here in sjsjwe employ a first order finite difference method to discretize and 
solve the integral equation ( |39| |. The solution to this equation contains high-order cancellations which can make it numerically 
challenging to solve. We have found that for a sharp fc-space filter this method is robust for all S and for arbitrary barriers 
providing that care is taken in the numerical evaluation. However, we find that the addition of perturbation terms can induce 
numerical inaccuracies that lead to fluctuation and divergence of the solution for 5 > 1. To circumvent this problem we 
proceed as follows. 

We begin by rewriting eqn. (411 in a more general form, considering all perturbation terms: 

AS 



' erf 



f B{S,) 



1 - erf 



fiSj) 



B{S,) - B{S,) 

V^iSj - SO 



-2Y,^''{Sj,S,))aS. 



(Bl) 



We now adjust the weight given to the flrst term on the right-hand side of eqn. (|41[) from AS/2 to AS /a. Solving for 



.f{Sj 



fiS, 



erf 



BjS,) 



j-1 



B{S,)-B{S.] 



2J2'^''{S,,S.) AS 



a 

AS' 



(B2) 



Note that f{Sj) is proportional to a. However, all /(5'i<j) also depend on a. In particular, if we make a slightly smaller 
than 2 we cause successive f(Sj) to oscillate around the true solution, effectively damping the numerical noise which otherwise 
causes divergence. 

Figures [B 1 1 and[B2l illustrate the convergence achieved using this numerical recipe. In both cases we compare our numerical 
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Figure Bl. The relative error of our numerical method with respect to the analytic solution for a sharp fc-space filter, Gaussian-random 
density field, and constant barrier. Results are shown for different values of the numerical parameter, a. The mesh size is fixed at 600 in 
all cases. 



results with the analytic solution for a sharp fc-space filter, a Gaussian random density field and a constant barrier. Figure [Bl] 
shows the relative error in the numerical and analytic solution for different values of a and a fixed mesh size of 600. While 
the error is lowest for a — 2 over most of the range, it is clearly growing for the largest values of S. We find that it continues 
to do so as S is increased further. As our solution is typically valid only for small S this may not be a problem. Nevertheless, 
we find that adopting a value of a slightly less than 2 results in acceptably small errors, and an error term that decreases 
with increasing 5*. Figure B2 considers different mesh sizes for fixed a = 1.5. The numerical solution is stable (with increasing 



accuracy) as the mesh size is increased, demonstrating that our solutions are robust and numerically converged. 
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Figure B2. The relative error of our numerical method with respect to the analytic solution for a sharp fc-space filter, Gaussian-random 
density field, and constant barrier. Results are shown for different mesh sizes (as indicated in the figure; here n is the number of mesh 
points in the range S = [0, 10]), and for a fixed a = 1.5. 



